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We analyze the configuration of a tliin rotating accretion disk, which is embedded in a magnetic 
field inducing a backreaction in the gravitating plasma. The aim of this study is to determine 
the conditions under which the gaseous accretion model of Shakura can be reconciled with the 
magneto-hydrodynamical picture requested to trigger the underlying turbulent behavior. We focus 
our attention to the generalized Ohm equation in order to understand if the plasma backreaction 
is able to provide the proper toroidal current, allowing a non-zero infalling velocity. In the limit 
of linear plasma backreaction, this analysis shows how the Shakura profile of accretion turns out 
to be inconsistent. In particular, comparing the azimuthal and the generalized Ohm equilibrium 
equations, we argue that it is not possible to maintain a constant rate of accretion. A non-stationary 
scenario for the disk configuration is then outlined and it results into a transient process which is 
however associated to a vanishing accretion rate. 
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INTRODUCTION 



One of the fundamental processes, governing the be- 
havior of a wide class of astrophysical sources, is the ac- 
cretion of matter onto a compact object The mate- 
rial surrounding massive accretors is, in general, in the 
state of a completely ionized plasma, as a consequence of 
both the high temperature and low density conditions. 
Nonetheless, the first proposal to describe the accretion 
process, due to Shakura in the early seventies 0,01; was 
based on a hydrodynamical approach for the infalling of 
a gas-like matter. This simplification was necessary to 
shed light on the mechanism underlying the angular mo- 
mentum transport. This aim was successfully reached 
identifying the shear viscosity, associated to the differen- 
tial rotation of the disk, as the driving phenomenon. 

The Shakura prescription is still widely adopted as a 
basic estimation for the accretion rate of astrophysical 
system versus viscosity. Indeed, despite recent multidi- 
mensional simulations ^-'t'I do not directly rely on the 
Shakura idea, nonetheless this model has still a timely 
profile both when theoretical Q and numerical [13] 
analyses are addressed. Since stellar accretion disks are 
mainly characterized by thin equilibrium configurations, 
the use of the one-dimensional model appears a well- 
grounded approximation and the comparison with de- 
tailed simulations is certainly viable, especially in view 
of the resulting values for the Prandtl number. The im- 
plementation of the Shakura prescription relies on its 
capability to offer a reliable paradigm for observations. 
Nonetheless, the microscopic features of the accreting 
plasma are typically compatible with a quasi-ideal sys- 
tem, for which the shear viscosity coefficient takes values 
much smaller than those ones requested to account for 



the observations. Shakura, being aware of this discrep- 
ancy, proposed that the dissipative effects arise from a 
quasi-ideal turbulent behavior, restated as an effective 
viscous laminar fiow. 

The solution to the puzzling picture concerning the 
linear stability of an axisymmetric rotating disk (incom- 
patible with the onset of turbulent regimes) was offered 
by the astrophysical implementation that Chaiidrasekar 
gave [ll| to the so-called magneto-rotational instability 
(MRI), just derived by Velikov ^ but in a slightly dif- 
ferent context. The merit of having understood how the 
MRI could be the proper mechanism to generate a dis- 
sipative dynamical regime, is due to Balbus and Hawley 
[13| (the analysis of alternative mechanisms responsible 
for unstable plasma modes in the disk can be found in 
Indeed, they clarified the role that an arbitrar- 
ily small magnetic field can play in developing an efficient 
turbulent behavior when coupled to a differentially rotat- 
ing plasma. The crossmatching of the Shakura idea of an- 
gular momentum transport with the MRI (here named 
the Standard Model for accretion), triggering of a tur- 
bulent regime, appears as a solid theoretical framework 
consistent with a significant set of observations [l|, [ll] . 

The present analysis is focused on a relic puzzle con- 
tained in such a picture (th e so-called anomalous resis- 
tivity of the disk plasma '1^) and on the suggestion that 
a revision of the basic hydrodynamical Shakura idea is 
mandatory. Indeed, including a magnetic field in the 
equilibrium configuration is a natural perspective when 
dealing with a real astrophysical source, often endowed 
with a significant magnetosphere, but it implies that the 
plasma (or simply the magneto-active fiuid) backreaction 
must be accounted too. In particular, due to the axial 
symmetry, the generalized Ohm law (GOL) requires that 
a relevant resistivity coefficient is postulated in order to 
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deal with a non-zero radial infall velocity. In fact, the 
current density arising in the plasma in view of its backre- 
action to the central object field (which has no associated 
intrinsic current since it is a vacuum field and it is approx- 
imated here by a dipole profile) can not be, in general, 
sufficiently strong to guarantee the balance accounting 
for the very small microscopic value of the plasma resis- 
tivity. In this paper, rather than investigate the mech- 
anisms generating the anomalous values of the resistiv- 
ity, we show that this puzzle enclose an inconsistency 
of the Shakura idea of accretion, when the full magneto- 
hydrodynamical (MHD) scheme is addressed. As natural 
development of our configuration scheme, we analyze a 
non-stationary profile of the disk which however turns 
out to provide a vanishing accretion rate on the compact 
object. 

The new feature of the one-dimensional model we 
present here (concerning also higher-dimensional numer- 
ical simulations) is the proof that the Shakura idea for 
the angular momentum transport holds only in the pres- 
ence of anomalous resistivity. This way, such a quan- 
tity becomes a theoretical prescription more than a phe- 
nomenological request. We reach this issue by a detailed 
balance of the magneto-static plasma backreaction and, 
apart from non-trivial technicalities, it is expected to 
hold in more than one-dimension too. The difficulty we 
outline here in reconciling the plasma magnetic backreac- 
tion with the Shakura prescription is particularly inter- 
esting in view of recent experimental 18] and numerical 
[igj issues which suggest the non-linear instability of neu- 
tral fluid in Keplerian rotation (see also [lOl); ^ '^^ll 
in view of the theoretical analysis which states (on 
a kinetic level) the stability of a magnetized plasma hav- 
ing a certain temperature anisotropy (against the MRI 
paradigm) . 

The necessity for a reformulation of the accretion pro- 
cess was proposed in [H,!!!] (see also [13]), where a paral- 
lelism between the axial symmetric configuration of Toka- 
mak fusion reactors and accretion disks was proposed. As 
underlying picture of the new point of view, in [25l [26j 
(see also the global extension provided in jlTf), it has 
been determined the existence of a crystal-like (i.e., a ra- 
dially oscillating profile) configuration for the magnetic 
field of a thin purely rotating disk. This scenario con- 
sists of a radial periodic structure of the magnetic flux 
surfaces, that, in the limit of a strong backreaction of the 
plasma, leads to a fragmentation of the profile into a ring 
series (for an analysis of jet formation in this scheme, see 
also [27[). The main merit of this new proposal is to 
outline how the Lorentz force, associated to the plasma 
backreaction, can play a relevant role in fixing the con- 
figuration, in view of the small-scale structure it acquires 
in the disk. More specifically, the Lorentz force is the 
proper contribution able to balance (already in the case of 
a linear backreaction) the higher order correction to the 
Keplerian centripetal force, in place of the gravitational 
effects. Indeed, the balance between the centripetal and 
gravitational forces takes place when identifying the an- 



gular velocity of the disk with its Keplerian profile at 
all orders of approximation in the mass density. There- 
fore, the additional term induced in the centripetal force 
by the backreaction must be included in the magneto- 
static equilibrium. Since also the Lorentz force is due 
to the plasma backreaction, it results natural to require 
that both the centripetal correction and the Lorentz term 
stand at the same equilibrium scale, as in ,25j . It is just 
this feature that makes our equilibrium configuration dif- 
ferent from previous solutions of the thin disk problem 
(see [H, [H liMIl). This is the line of thinking at the 
base of the present treatment for accretion process. 

The paper is organized as follows: in Section [TTl we dis- 
cuss the fundamental framework in which our analysis is 
performed. We introduce the fundamental quantities and 
the concept of vertical average necessary to reduce the 
dynamics to a one-dimensional problem. In Section IIIIl 
we construct the fundamental equations governing the 
equilibrium configuration. In particular, we recognize in 
the radial equation a separation between the hydrody- 
namical and the magnetic backreaction scales. Section 
llVl is devoted to outline the inconsistency of the Shakura 
model in the presence of a magnetic field by a detailed 
analysis of the azimuthal momentum conservation equa- 
tion. We also show how an analytical solution, radially 
oscillating and damped in time, can be obtained but it 
corresponds to a net zero accretion rate of the disk. We 
then conclude this section providing an estimate showing 
how the only viable approach to the Shakura prescription 
implies an anomalous resistivity coefficient. In SectionlVl 
we present a physical discussion concerning the nature of 
the magnetic backreaction at the ground of our approach. 
An astrophysical characterization of our model is offered 
in order to focus the nature of the obtained inconsistency. 

II. THE UNDERLYING PARADIGM 

We now construct a non-stationary MHD model for 
the equilibrium of a thin accretion disk, surrounding a 
compact astrophysical object. In axial symmetry, the 
magnetic field B reads as 

B = -Brdy^ip/r + e,^drik/r , (1) 

in view of its divergenceless nature ((r, 0, z) being cylin- 
drical coordinates). Here, we set to zero the toroidal 
component i?^, while ip — 'ip{t,r,z'^) denotes the mag- 
netic flux surface. Let us now split this function as 

^{t,r,z^)~Mr,z^)+Mr,t,z^) , (2) 

denoting by V'o the flux surface associated to the back- 
ground magnetic field (proper of the central object), and 
by ipi the contribution due to the backreaction of the 
disk plasma. In what follows, the flux surface of the 
background field is taken in a reliable dipole form and 
we assume a linear response of the plasma ipi ipo, for 
which the induced magnetic field remains small with re- 
spect to the background one, i.e., drtpo ^ dripi- The 
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density current can be split in the form J = Jo + Ji, 
where Jq = 0, as requested by the vacuum nature of 
the corresponding field. Furthermore, the z-dependence 
of tpi can be neglected when calculating the backreac- 
tion component Ji of the current. In fact, considering 
the regime where the /3-parameter of the plasma is much 
greater than unity, namely /? ~ /u^ 3> 1 (ws and va be- 
ing the plasma sound and Alfven velocity, respectively), 
the typical vertical scale of the configuration 77 ~ 
(where H — H{r) is the half-depth of the disk) is much 
greater than the corresponding radial one A ~ va/^ [111 
(the consistency of this scheme can be verified a posteriori 
from Eqs.([2T|)). This hierarchy of the typical configura- 
tion scales is naturally expected to be maintained in the 
linear backreaction regime and it allows to neglect the 
contribution that the z-derivatives of tpi would provide 
to the toroidal current density. 

In agreement to the Standard Model for accretion [l|, 
let us now consider an essentially one-dimensional ax- 
isymmetric scheme, in which the z-dependence is aver- 
aged out using 

^'^'^'*^"2]/ j7i^^t^^^)dz, (3) 

where / denotes a generic coordinate function and we dis- 
regarded odd z-quantities because their average vanish 
by definition. This simplification is due to the thinness 
of the disk and the z-derivatives of all the even quanti- 
ties with respect to the equatorial plane can be neglected 
together with the vertical component Vz of the fluid veloc- 
ity (vanishing on the equatorial plane) . Accordingly, the 
mass density p can be replaced by the superflcial density 
defined by T,(t,r) = 2Hp, where the thinness condition 
H/r <^ 1 holds. While, averaging out the configuration 
equations, we replace the radial velocity Vr by a mean 
value Ur weighted over the density and the same for its 
derivatives (or products of these quantities), i.e., 

_ 1 1 f^" 

= V / P^rdz , dt.r Ur = — p dt.r Vrdz . (4) 

respectively. Another important property characterizing 
the disk configuration is its angular velocity $7, taken 
equal everywhere in the disk to its equatorial value. 

Introducing the disk accretion rate M = —2TrrT,Ur [l|, 
the continuity equation writes 

QtS - drM/2TTr = . (5) 

We note that, while d^Vz does not vanishes on the equato- 
rial plane, the total z-derivative in this equation provides 
a very small edge term after integration. 

After averaging over the vertical profile, the only sur- 
viving component of the magnetic field B is 

B = ezdri^/r, i/i(t, r) = V^o -f V-i , (6) 

where tjjo = n/r (with p, = const.), coherently with the 
dipole expression. Averaging now the current density 



J = Ji, we stress that the key point of the proposed 
scheme consists in the large value that the resulting back- 
reaction current component, i.e., 

— c/1— \ c — 

= dr[- dr^i ) dli,x , (7) 

takes in view of the small scale A of the plasma backre- 
action (in fact, in the last passage above, we made use of 
the inequality A -C r so neglecting the derivative of the 
l/r factor in front of dr'ipi). 

III. FUNDAMENTAL EQUATIONS 

Let us now consider the system provided by the Fara- 
day equation and GOL, i.e., 

dtB + cV X E = , (8) 
ca E + au X B - cJ ^0 , (9) 

respectively [13 , where E denotes the electric field and a 
the electric conductivity. The system above can be easily 
combined using the expression ^ for the magnetic field. 
In particular, the z-component of Eq.® provides Ez = 0, 
while the radial one writes 

Er = -dr<P = -i^/c) dri) , (10a) 

where $ denotes the electrostatic potential. In this 
scheme, the z-component of Eq.® yields to the relation 
E^ = —dtipi/cr, while the two other ones identically van- 
ish. Finally, the toroidal component of Eq.® writes 

9t?/5i -|-Ur9rV'o = -cr Ji(0)/cr , (10b) 

We stress that Eq. ljlObp is equivalent, apart from a global 
radial derivative, to the z-component (the only non- 
vanishing one) of the magnetic-field induction equation 
obtained substituting Eq.® into Eq.®. 

We observe that Eq. ljlOap is compatible with the as- 
sumption = Q{ip(t,r)). Since in the absence of back- 
reaction the disk rotation is Keplerian, i.e., 

no = nK = VGMs/r^ = ^Jgms^I/p^ , (11) 

[Ms being the mass of the central object), the simplest 
expression satisfying the corotation condition above and 
able to recover the Keplerian nature at the lowest order 
can be postulated as 

n{^) = ^JgMs^Vp^ ^nK(l + 3^i/2Vio) , (12) 

where it is essential that the amplitude of the perturbed 
magnetic surfaces remains always smaller than the back- 
ground value ijjQ — p/r'^. In the same approximation 
scheme, setting rjB — (? j^Tia and using the relation be- 
tween Ur and M, Eq. ljlObp rewrites 

E^tV^i + Mp/2'Kr^ = ErysS^V'i • (13) 
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The radial equilibrium equation [l| can be split into 
three components, corresponding to different orders of 
the equilibrium configuration. The first two naturally en- 
ter the Shakura model for a non-magnetized fluid, while 
the third one is due to the plasma backreaction and it 
balances the corrections of the centripetal force by the 
Lorentz one. They write 



no 

Ti{dtUr + UrdrUr) 



VtK , (14a) 
-dr^ , (14b) 
-{3GM,/^^H)i:r^rPi , (14c) 



respectively, where 11 = 2Hp {p denoting the thermo- 
static pressure). The possibility to restate the radial 
equilibrium in these three separated force balance equa- 
tions is a fundamental step of our approach. While the 
Keplerian nature of the background configuration is a 
well posed assumption for a thin disk, see [321 > the split- 
ting between the hydrodynamical term and the mag- 
netic backreaction can take place only if the plasma /3- 
parameter satisfy the condition (3 <^ L/\ (hence X <^ L) 
being L the variation scale of the pressure. Furthermore, 
using M and Eq.dU, Eq. (|14b|) rewrites 



dtM + dr(Mur) = 2TrrdrIl 



(14d) 



The vertical equilibrium equation [l[ fixes, via the 
gravostatic equilibrium, the z-dependence of the mass 
density p. The features of the vertical equilibrium, sen- 
sitive to the equation of state, are not needed for the 
addressed task, apart from fixing H{r) ~ Vs/^k ^ r. 

Finally, using the expression (|12p for 51, the azimuthal 
equilibrium equation [l| writes, in the presence of a mag- 
netic field, 



(15) 



where rjy denotes the shear viscosity coefficient associated 
to the disk differential rotation (we have neglected the 
gradient of ipi with respect to its second r-derivatives). 



IV. INCONSISTENCY OF THE MAGNETIZED 
SHAKURA MODEL 

We remind that the Shakura model [3] is based on the 
idea that the accretion rate of the disk is driven by the 
background shear viscosity. Requiring that the accretion 
rate is provided by the usual Shakura expression, i.e., 
Msh — GT^Hrjy [l|, indeed we are implementing the idea 
that the magnetic field of the central object has the only 
role to allow for the MRI. This way, the Shakura model 
is preserved because the turbulence arising from the MRI 
induces the effective values of the dissipation coefficients 
on which the phenomenological predictability of the ac- 
cretion rate Msh is based. The source of inconsistency 
can be easily identified in the new equation arising in the 



MHD scenario, i.e., GOL, especially its azimuthal com- 
ponent. It is just this new equation which can not be 
reconciled with the azimuthal component of the momen- 
tum conservation (jlSp . once the Shakura accretion rate 
is imposed. 

In fact, following this prescription in order to recover 
the standard behavior, we assume the decomposition of 
Ea. (|15p in two equations, one for the background regime 
and the other one for the backreaction response, i.e., 

Mdrir^nK)/2T: + dr {2Hr]yr^drnK) = , (16a) 
E dtipi - 2Hr)y dl^i = , (16b) 

respectively. Eq. (|16ap matches the standard Shakura be- 
havior for the angular momentum transport. In fact, if 
we focus on the stationary configuration by setting to 
zero all the partial time derivatives, it reduces to the ordi- 
nary Shakura relation Msh = QTrHrj^. Moreover, Eq.® 
gives the fundamental constraint M = const. 

In the stationary regime, the inconsistency of the Stan- 
dard Model for accretion disks comes out immediately: 
Ea. (ll6bp requires d'^ipi = implying, by virtue of GOL 
([T3|. M ~ 0, i.e., a vanishing disk accretion rate. By 
other words, in the presence of a magnetic field, a back- 
reaction current of the plasma is mandatory in order to 
balance the accretion rate in Eq. (jl3p . but such a current 
must vanish (or, more realistically, it must be negligible) 
to safe the accretion-rate constant profile of the Shakura 
prescription. This type of inconsistency can not be re- 
moved even retaining the z-dependence of tpi . 

We now focus on the non-stationary sector inferring 
that the background term proportional to viscosity can 
be neglected, while the accretion process should be driven 
by the toroidal density current. Thus, Eg. dTSI) rewrites 

E^tV-i - M^^/6^Tr^ ^ 2Hr]yd^'iPi , (17) 
and combining Eas. (fT3)) . (|14cP and (IT7)) . we get 

2Mp./3nr^ = (2Hr]y - j:r]B)SGMsT.r^ ipi/Hp.'^ , (18) 
dti^i = -{2Ht], + j:rjB/3)9GMsr'' i^i/4H^i^ . (19) 

The dynamical system involves four unknowns E, M, 
-01 and n obeying Eqs.©, (|14dl) . pSj) and (HH) (com- 
pleted retaining Eq. (|14cP too). In the case of linear 
backreaction, the superficial density must admit the per- 
turbed decomposition E = Eo(?') -f- Ei(t,r) (Ei ^ Eq) 
where Eg is associated to the background profile, while 
El is the small time-dependent correction due to the 
plasma backreaction. In order to get a physical insight on 
the dynamical content of the accretion paradigm traced 
above, we can make some reasonable simplifying assump- 
tions setting 
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E = Eo = 2Hm/r 
rjB — const. , 
2Hriv = Y.r]B/np , 



(20a) 
(20b) 
(20c) 
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where m = const, and np = const denotes the Prandtl 
number. The behavior of rpi is now described by 
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(21a) 
(21b) 



= 6GMsm/fi^ . (21c) 

The solution reads as 

^1 = ^1 sin(fcr) exp [ - (i + ^) yy^fc^ t] , (22) 

where ipi is an integration constant and substituting this 
expression into Eg. dTSl) . we finahy get 



M = {l/np - 1) {3TT/fi)mHriBk'^ ipi 



(23) 



Although the resulting M is comparable in amplitude to 
Msh, we stress that the radial infalling velocity is de- 
scribed by the behavior Ur oc —r'^'ipi. Since kr ^ 1 
by hypothesis {k ~ 1/A), the radially oscillating term 
of ipi varies much rapidly than and, therefore, the 
last can be approximated as constant when estimating 
the accretion profile. As a consequence, the obtained 
non-stationary model provides a radial infalling velocity 
having a vary fast radial sinusoidal behavior damped in 
time by the resistive effects. Clearly, such a profile is not 
associated with a real infalling of material on the cen- 
tral object, but it is essentially a local accretion process 
which accumulates material both from outside and inside, 
around a fixed radial coordinate. The same scenario has 
been also outlined in [s^ [sj] within a two-dimensional 
local configuration. 

Our work claims the impossibility to reconcile the 
Shakura Model with a fully consistent MHD scheme in 
which the plasma backreaction can not be neglected. The 
open question is then if the assumption kr ^ 1 {i.e., 
A ^ r) at the base of our analysis is really obliged or 
it can be somehow avoided. We elucidate this point by 
a chain of estimates based on the fundamental equilib- 
rium equations. In the stationary regime, if we adopt the 
Shakura relation for the viscosity r]^ = 2aTiQVs/3 (where 
a < 1 is a parameter of the model), Msh provides the 
radial infalling velocity Ur ^ aVgH/r. At the same time, 
Eq.® gives Ur ^ Bic^ / B^aX where Bq and Bi are the 
background and backreaction magnetic field amplitudes, 
respectively, and we recall that A denotes the magnetic 
backreaction scale. Using now the definitions of a and 
•qB , and comparing the obtained expressions of the radial 
velocity, one can easily get 



Bi r 
B^ A 
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(24) 



Thus, in the linear regime (Bi <C Bq), if we require that 
the backreaction currents induced in the plasma are not 
relevant, i.e., A ~ r, in order to avoid a vanishing accre- 
tion rate in the Standard Model, one must address huge 
values of the resistivity coefficient since rip > 1. 



This work demonstrates that the Shakura model is in- 
consistent in the presence of a small-scale magnetic back- 
reaction able to account for an enough strong toroidal 
current. Indeed, only such a current allows a balance 
of the electron force in the presence of quasi-ideal (ki- 
netic) values of the resistivity coefficient. If we try to 
avoid this inconsistency by considering large-scale back- 
reaction, then we have to deal with a very small toroidal 
current which can be neglected everywhere in setting the 
equilibrium configuration, but in the toroidal GOL. In 
fact, in such an equation the term t/b Ji(0) must be nec- 
essarily retained to balance the Lorentz force acting on 
the electrons. As a consequence, if the backreaction is 
of very large-scale and the corresponding toroidal cur- 
rent are tiny, the only way to preserve the consistence 
of the stationary GOL is to consider a huge value of the 
coefficient r]B, i.e., to deal with the so-called anomalous 
resistivity. This statement is well-elucidate by the esti- 
mate above. 

Putting together the present analytical treatment, 
which leads to the solution (22) , with the qualitative fea- 
tures outlined by the estimate (|24| . we are able to get 
light on the origin of the anomalous resistivity puzzle: 
when the numerical simulations are compared with the 
observations, a strong value of the resistivity coefficient 
is needed because very small-scale magnetic backreaction 
is forbidden by the compatibility of the equilibrium con- 
figuration. 



V. PHYSICAL HINTS AND FINAL REMARKS 

It is now worth clarifying the real nature of the plasma 
backreaction we are referring to and which is the phys- 
ical meaning of the solution (f^. Indeed, the presence 
of the external magnetic field is responsible for a series 
of correlated non-linear processes whose net effect is the 
emergence of the effective visco- resistive properties ^2^ . 
Such phenomena are called plasma backreaction too and 
it is crucial to stress how it intrinsically differs from the 
disk reaction we are accounting for. 

The fundamental feature to be focused in the plasma 
backreaction is the interplay which takes place in the 
disk between the MRI-triggered turbulence (responsible 
for the magnetic field enhancement) and the turbulent re- 
connecting transport (transferring outward such a field) 
[35l - l38| . By other words, the dynamo effect produces 
a growth of the magnetic field amplitude, which would 
stop the MRI instability, but the latter is continuously 
restarted by the decrement of the magnetic field due to 
the turbulent transport of its fiux in the outer regions of 
the disk, and in the disk corona. We can estimate the 
typical time scale of these processes in terms of the char- 
acteristic growth rate for the MRI instability which, in 
the best case, is tmri ~ ^/^k IM- Finally, it is worth 
noting that the essential background for such instabil- 
ity is the purely rotating Keplerian disk, lacking of any 
accretion feature. 
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The backreaction scheme we present here, and which 
is responsible for the inconsistency of the Shakura pre- 
scription in the MHD sector, has a very different ori- 
gin. In fact, it takes place on much shorter time scale 
than the MRI triggered turbulence and it arises nearby 
a background configuration which, in addition to the 
stripped rotating disk possesses, possesses an hydrody- 
namical scale accounting, in principle, for an accretion 
process. By other words, the net effect of the compli- 
cated interplay between the turbulent dynamo and re- 
connecting transport is just the emergence of the non- 
zero accretion of the disk, that we assume at the ground 
of our analysis. In this sense, the above mentioned pro- 
cesses are here summarized in the non-zero viscosity and 
resistivity coefficients. 

The characteristic time scale of the backreaction we 
analyzed, as it results from equation (|22|) . is of the order 
tbr ~ np/r]Bk^, where we excluded the limit np ^ 1. 
Since we are implementing our backreaction on a stan- 
dard (effective) visco-resistive scenario, we can adopt the 
value of rjB when the Shakura prescription for the vis- 
cosity coefficient is considered, i.e., r]B ~ npaVgH ~ 
npaH-flx ■ Hence, we get the fundamental relation 
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I 



tbr 



auJK{kHY a{kH)' 



tmri 



(25) 



Although the parameter a can take small values and 
its range of variation is from 0(10"'^) to 0(1), the quan- 
tity kH can take huge values in the disk. An estimation 
for the wave-number (|21c|) has been provided in [13] and 
it gives a minimum value of about k,nin ~ lOifm"^. 
Thus, if we take for the disk half-depth the typical value 
H ~ 10'^ Km, even in the best case a ~ 10~^, we get 
the relation tbr ^ 10~^tmri ■ Such a striking difference 
in the characteristic time scale of the MRI triggered tur- 
bulence and the backreaction here analyzed outlines the 
correspondingly different nature between these dynami- 
cal processes. The very extremely small period of time, 
over which the backreaction driven by ■01 is damped, 
confirms the inconsistency of the Shakura model in the 
presence of a magnetic field. In fact, there is no way 
to deal with a longstanding current in the disk, able to 
account for the proper balance of GOL. This current is 
induced by the plasma disk backreaction and it is al- 
most instantaneously damped by the dissipative effects 
(the net result of the turbulent processes on which the 



Shakura paradigm relies on). As a consequence of this 
picture, the solution (|22p more than representing a phys- 
ical regime of the disk dynamics (to be analyze in view 
of its possible implications), it results a rapid damping 
profile of any radial oscillating behavior arising in the 
disk plasma by means of the balance between the in- 
duced centripetal and Lorentz force, respectively. In this 
sense, the crystalline morphology of the disk predicted 
in [25[ is just a spatially short-scale and very rapid tran- 
sient phenomenon. This is the base of our criticism to the 
Shakura prescription, because if the magnetic backreac- 
tion is damped by the dissipation required to deal with 
a satisfactory accretion, where does the current able to 
ensure the proper azimuthal GOL come from? 

We conclude observing the the solution to this puzzling 
picture can not come from the interaction of the disk 
with its environment, for instance with the disk corona 
[39I l40j . In fact, the process responsible for exchange 
of energy between the disk and its corona seems to be 
driven, on a magnetic level, on very large spatial scales 
411 . The backreaction we studied here, in the extreme 
cases, is a phenomenon having the spatial scale of the me- 
ters and the typical duration of /iseconds. It is just this 
peculiar feature which allows us to claim firmly the pres- 
ence of an inconsistent scheme, as soon as, the hydrody- 
namical and magnetic scale are separated to recover the 
Shakura profile. Indeed, it is clear that such a scheme 
can not be affected by the typical space and time scales 
of the fundamental physical processes governing the disk 
profile evolution. Therefore, only two possible scenarios 
remain viable in the perspective of the present issue: (i) 
it requires an enormous anomalous value of the resistivity 
coefficient, as discussed above; (ii) the separation of the 
hydrodynamical and the magnetic scales is removed and 
a deep interchange between the kinetic and magnetic en- 
ergy of the disk takes place (the disk morphology is then 
rather far from the Shakura prescription). This second 
perspective requires the set-up of sophisticated numeri- 
cal simulations, able to access the very short space and 
time scale, whose contribution can not be, in principle, 
disregarded as outlined here. 
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